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Abstract 

Numerous exact relations exist that relate the effective elastic properties of composites to the elastic properties of their com- 
ponents. These relations can not only be used to determine the properties of certain composites, but also provide checks 
on the accuracy on numerical techniques for the calculation of effective properties. In this work, some exact relations are 
discussed and estimates from finite element calculations, the generalized method of cells and the recursive cell method are 
compared with estimates from the exact relations. Comparisons with effective properties predicted using exact relations show 
that the best estimates are obtained from the finite element calculations while the moduli are overestimated by the recursive 
cell method and underestimated by the generalized method of cells. However, not all exact relations can be used to make such 
a distinction. 

1 Introduction 

Exact relations for the effective elastic properties of two-component composites can be classified into three types. 
The first type consists of relations that have been determined from the similarity of the two-dimensional stress and 
strain fields for certain types of materials. These exact relations are called duality relations HI. The second type of 
exact relations, called translation-based relations, state that if a constant quantity is added to the elastic moduli of 
the component materials then the effective elastic moduli are also "translated" by the same amount. Microstructure 
independent exact relations, valid for special combinations of the elastic properties of the components, form the 
third category [2J. The known exact relations are directly applicable only to a limited range of properties of the 
components. Therefore the utility of these relations lies not only in determining the effective elastic properties of 
a small range of composites but also in evaluating the accuracy of numerical and analytical methods of computing 
effective properties. In this work, predictions from exact relations are compared with estimates from finite element 
calculations, the generalized method of cells (GMC) [3], and the recursive cell method (RCM) The goal is 
to assess the effectiveness of these relations in evaluating the accuracy of the three numerical methods, especially 
with regard to high modulus contrast materials such as polymer bonded explosives. 

Five exact relations are explored in this work. The first is a duality-based identity for the effective shear 
modulus that is valid for phase-interchangeable materials [5 ]. The second is a set of duality relations that are valid 
for materials that are rigid with respect to shear [1J. Two translation-based relations are explored next - the CLM 
theorem [6] and a relation for symmetric composites with equal bulk modulus 13. The microstructure independent 
Hill's relation Q is explored last. 

2 Phase interchange identity 

A symmetric composite is one that is invariant with respect to interchange of the components. A checkerboard, 
as shown in Figure [TJ is an example of a symmetric composite. The phase interchange identity Q for the effec- 
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Figure 1 : Representative volume element for a checkerboard composite. 

tive shear modulus of a symmetric two-dimensional two-component isotropic composite is a duality-based exact 
relation that states that 

Geff = ^GTGs (1) 

where G\, G2 are the shear moduli of the two components and G e ff is the effective shear modulus. 

The phase interchange identity is valid only for isotropic composites. In a finite-sized representative volume 
element (RVE) for a checkerboard composite the shear modulus is not the same all directions and hence isotropy is 
not achieved. The two-dimensional stress-strain relation for such a RVE with "square symmetry" can be written as 
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where an, 022, 012 are the stresses; en, 622, £12 are the strains; K is the two-dimensional bulk modulus, ji is the 
shear modulus when shear is applied along the diagonals of the RVE, and ju 2 is the shear modulus for shear along 
the edges of the RVE. 

The numerical verification of the phase interchange identity, therefore, requires that the components of the 
composite be chosen so that the difference between [i 1 and fi 2 for the composite is minimal. This implies that the 
components should have a weak modulus contrast. 

Numerical estimates of the effective elastic properties of the checkerboard composite shown in Figure [T] were 
obtained using finite elements (FEM), the recursive cell method (RCM) and the generalized method of cells (GMC). 
Following the requirement of low modulus contrast, both components were assigned a Young's modulus of 15,300 
MPa. The Poisson's ratio of the first component was fixed at 0.32 while that of the second component was varied 
from 0.1 to 0.49. The FEM calculations were performed using a mesh of 256x256 four-noded square elements. 
The RCM calculations used a grid of 64x64 subcells with blocks of 2x2 subcells and each subcell was modeled 
using one nine-noded element. The GMC calculations used 64 x 64 square subcells to discretize the RVE. 

Figure [2] shows a comparison of the exact effective shear modulus for the checkerboard composite with esti- 
mates of /x 1 and /x 2 from the three numerical approaches. The results show that all the three methods perform well 
(the maximum error is 0.1%) in predicting the effective shear modulus when the modulus contrast is small, i.e., 
when the composite is nearly isotropic. It can also be observed that the values of /j, 1 and /x 2 are within 1% of each 
other for the chosen component moduli. 
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Figure 2: Validation of FEM, RCM and GMC using the phase interchange identity for a checkerboard composite. 

2.1 Range of applicability 

The question that arises at this point is whether the three numerical approaches can predict the phase interchange 
identity for larger modulus contrasts. Numerical calculations have been performed on the checkerboard microstruc- 
ture to explore this issue. The first component of the checkerboard was assigned a Young's modulus of 15,300 MPa 
and a Poisson's ratio of 0.32. For the second component, the Poisson's ratio was fixed at 0.49 and the Young's 
modulus was varied from 0.7 MPa to 7000 MPa. 

Figure |3] shows plots of the effective fi 1 and fj, 2 versus shear modulus contrast for a checkerboard RVE. The 
plots confirm that when the modulus contrast between the components of the checkerboard exceeds 2, the material 
can no longer be considered isotropic since the values of fi 1 and fj 2 are considerably different from each other. 
However, the values of /j, 1 predicted by FEM are quite close to the effective shear modulus G e ff predicted by the 
phase interchange identity. This result suggests that the simulation of a diagonal shear may not be necessary to 
predict the effective shear modulus of an isotropic composite when the finite element approach is used. It also 
implies that the phase interchange identity can be used for a much larger range of modulus contrasts. The effective 
shear moduli predicted by GMC are considerably lower than that from the exact relation while the values from 
RCM are consistently higher. The RCM estimates worsen with increasing modulus contrast. If only the value of 
/1 1 is examined, the phase interchange identity indicates that the FEM approach is much more accurate than the 
GMC and RCM approaches. However, it is difficult to choose between GMC and RCM for high modulus contrasts 
composites. While the exact value of [i 1 is 10 times the value predicted by GMC, the corresponding RCM estimate 
is 10 times the exact value. These results confirm the findings of detailed numerical studies on high modulus 
contrast, high volume fraction polymer bonded explosives [8j[9j|4]]. 

2.2 Convergence of FEM calculations 

The checkerboard material provides an extreme case to test the convergence of the FEM solution because the corner 
singularities lead to high stresses that can only be resolved with refined meshes. Figure [4] shows the convergence 
of the effective fj, 1 and fj 2 with increasing mesh refinement for a checkerboard with a shear modulus contrast of 
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Figure 3: Variation of effective shear moduli with modulus contrast for a checkerboard composite. 

about 25,000. The effective fi 1 converges to a steady value when 128 x 128 elements are used to discretize the RVE. 
The shear modulus fj, 2 reaches a steady value when 256x256 elements are used. This is why the finite element 
calculations in this work were performed using 256x256 elements or more. 

RCM uses a finite element approach to homogenize blocks of subcells. When blocks of 2x2 subcells are used, 
some of these blocks can resemble checkerboards - especially at the first level of recursion for a two-component 
composite. The finite element convergence result suggests that RCM may overestimate the effective shear moduli 
by a factor of two if a block of four subcells is simulated using only four finite elements. 

3 Materials rigid in shear 

The stress-strain response of two-dimensional materials that are rigid with respect to shear can be represented by 
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where an, 022, 012 are the stresses; en, £22 and 612 are the strains, and S%j are the components of the compliance 
matrix. 

Two duality-based exact relations that are valid for two-component composites composed of such materials 
are ffl: 

Relation RSI If SnS 22 — (£i2) 2 = A for each phase (where A is a constant), then the effective compliance tensor 
also satisfies the same relationship, i.e., SffSf^ — (Sf 2 ) 2 = A e ff. This relation is true for all microstructures. 

Relation RS2 If the compliance tensors of the two phases are of the form Si = ot\ A and S2 = cx 2 A. where A 
is a constant matrix, then the effective compliance tensor of a checkerboard of the two phases satisfies the 
relation detS eff = SfiSg - {Sf 2 f ) 2 = ai a 2 (A n A 22 - {A 12 ) 2 ). 
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Figure 4: Convergence of effective moduli predicted by FEM with increase in mesh refinement for a checkerboard 
composite with a shear modulus contrast of 25,000. 



3.1 Relation RSI 

Figure [5] shows a square array of disks occupying an area fraction of 0.7. Numerical experiments have been 
performed on this array of disks to check if Relation RS 1 can be reproduced by finite element analyses, GMC 
and RCM. The S matrices that have been used for the disks (superscript 1) and the matrix (superscript 2), and the 
corresponding values of A are shown below. These matrices have been chosen so that the value of A is constant. 



Si 



1000 -300 
-300 1000 
0.001 



A = 9.1 x 10 5 , 



and 



1094.3 
-536.21 




-536.21 
1094.3 
0.001 



A = 9.1 x 10 5 



The shear modulus for both materials is 1000 (arbitrary units) - around 10 6 times the Young's modulus. Higher 
values of shear modulus have been tested and found not to affect the effective stiffness matrix significantly. 

Table [l] shows the values of S'ff, Sff and A e g calculated using finite elements (350x350 elements), GMC 
(64 x 64 subcells) and RCM (256x256 subcells). The ratio of the calculated A e g to the original A are also shown 
in the table. The modulus contrast between the two components of the composite is small, so the calculated 
effective properties are expected to be accurate (based on the results on the phase interchange identity for shear 
moduli). However, the results in Table [T] show that all the three numerical methods predict values of A e ff that are 
around half the original A. These results imply that all three methods (FEM, GMC and RCM) overestimate the 
effective normal stiffness of the array of disks. Relation RSI for materials rigid in shear may therefore be a very 
sensitive test of the accuracy of numerical methods even though the modulus contrast that can be used is small. 
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3.2 Relation RS2 



The second duality relation for materials that are rigid in shear requires (Relation RS2) is valid for the checkerboard 
geometry shown in Figure [T] The following values of the elastic properties have been used to test the accuracy of 
FEM, RCM and GMC in predicting this relation. 
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The duality relation requires that the effective compliance matrix of the checkerboard composite should be such 
that 

det(Seff) = Sffsg - (Sit) 2 = 9-10 x 10 6 . 

The FEM calculations were performed using 350 x 350 four-noded elements, the RCM calculations used 64 x 64 
subcells (blocks of 2x2 subcells) and the GMC calculations used 64x64 subcells too. The results from these three 
methods are tabulated in Table [2] The finite element calculations lead to quite an accurate effective compliance 
matrix and the deviation from the exact result is only around 25%. The GMC calculations overestimate the com- 
pliance matrix and the determinant of the compliance matrix is around 2.3 times higher than the exact result. On 
the other hand, the RCM calculations predict a compliance matrix that has a determinant that is only around 20 



4 The CLM theorem 

The Cherkaev, Lurie and Milton (CLM) theorem is a well known "translation" based exact relation for two- 
component planar composites |6|. For a two-dimensional two-component isotropic composite, this theorem can be 
stated as follows. 

Let the isotropic bulk moduli of the components be K\ and Ki. Let the shear moduli of the two components be 
G\ and G2. The effective bulk and shear modulus of a two-dimensional composite made of these two components 
are K c q and G e ff respectively. Let us now create two new materials that are "translated" from the original compo- 
nent materials by a constant amount A. That is, let the bulk and shear moduli of the translated component materials 
be given by 

1/Kf = \jK x - A ; 1/Kj = 1/K 2 - A ; 
1/Gj = 1/Gi + A ; l/Gl = 1/G 2 + A . 

The CLM theorem states that the effective bulk and shear moduli of a two-dimensional composite of the two 
translated materials, having the same microstructure as the original composite, are given by 

= - A ; 1/G5r = 1/GW + A. (4) 



Table 1: Two-dimensional effective compliance matrix for a square array of disks. 





ceff 
D n 


ceff 
°12 


A eff (Xl0 5 ) 


Aeff/A 


FEM 


850.35 


-536.32 


4.35 


0.48 


RCM 


847.25 


-538.17 


4.28 


0.47 


GMC 


871.75 


-517.14 


4.93 


0.54 
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Table 2: Effective compliance matrix for a checkerboard composite with components rigid in shear. 





S$ 


ceff 
D 12 


det(S eff )(xl0 6 ) 


det(S ett )/ai«2 det A 


FEM 


3282 


-2004 


6.75 


0.74 


RCM 


1655 


-7090 


2.23 


0.24 


GMC 


5007 


-2146 


2.05 


2.25 



The requirement of isotropy can be satisfied approximately by choosing component material properties that are 
close to each other. Since our goal is to determine how well GMC and RCM perform for high modulus contrast, 
choosing materials with small modulus contrast is not adequate. Another alternative is to choose a RVE that 
represents a hexagonal packing of disks. However, such an RVE is necessarily rectangular and cannot be modeled 
using RCM in its current form. It should be noted that RCM can easily be modified to deal with elements that are 
not square and hence to model rectangular regions. 

Another problem in the application of the CLM theorem is that the value of A has to be small if the difference 
between the original and the translated moduli is large and vice versa. If the value of A is small, floating point 
errors can accumulate and exceed the value of A. On the other hand, if A is large, the original and the translated 
moduli are very close to each other and the difference between the two can be lost because of errors in precision. 
Hence, the numbers have to be chosen carefully keeping in mind the limits on the value of the Poisson's ratio. 

The translation relation has been tested on the square array of disks occupying a volume fraction of 0.70 from 
Figure[5j This RVE exhibits square symmetry, i.e., the shear moduli // and fi 2 shown in equation Q are not equal. 
A unique value of the effective shear modulus cannot be calculated for this RVE. Instead, he value of the effective 
translated shear modulus is calculated from equation Q by first setting G e ff equal to fi 1 and then to /i 2 . These 
"exact" values are compared with the fi 1 and fi 2 values predicted using finite element analyses, GMC and RCM. 

The original set of elastic moduli for the RVE is chosen to reflect the elastic moduli of the constituents of poly- 
mer bonded explosives. These moduli are then translated by a constant A = 0.001. The original and the translated 
constituent two-dimensional moduli are shown in Table [3] (phase 'p' represents the particles and phase 'b' repre- 
sents the binder). It can be observed that the translation process creates quite a large change in the bulk modulus of 
the particles. Table|4]shows the effective bulk and shear moduli of the original and the translated material calculated 
using finite elements (350x350 elements), GMC (64x64 subcells) and RCM (256x256 subcells). The values of 
A eiT shown in the table have been calculated using the equation 

A err = (A/0.001 - 1) x 100, 
A = - l/Kj n = 1/4P - 1//4. 

Even though the modulus contrast between the two components of the composite is high, the effective properties 
predicted by FEM, GMC and RCM are close to each other in magnitude. The effective moduli of the translated 
composite are also quite close to that of the original composite as predicted by the CLM condition. The interesting 
fact is that all the three methods satisfy the CLM condition and the error is small (as seen by the values of A err . Of 
the three methods, FEM and GMC produce the least error while RCM produces the most error. 

Table 3: Original and translated two-dimensional constituent moduli for checking the CLM condition. 





K p 


G p 


Kb 


Gb 


K p /K b 


Gp/Gb 




(xlO 2 ) 


(xlO 2 ) 






(xlO 2 ) 


(xlO 2 ) 


Original 


9.60 


4.80 


10.07 


0.20 


0.95 


23.8 


Translated 


240.0 


3.24 


10.17 


0.20 


23.5 


16.1 
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5 Composites with equal bulk modulus 



The translation procedure can also be used to generate an exact solution for the effective shear modulus of two- 
dimensional symmetric two-component composites with both components having the same bulk modulus [5 ]. This 
relation is 

K eff = K = K\ = K 2 

G eff = , K (5) 

-W( 1+ £>( 1+ £) 

This relation has been tested on the checkerboard model shown in Figure [T] using the component material properties 
given in Table [5] The exact effective properties for the composite, calculated using equation ([5]), are also given in 
the table. The values of the effective moduli calculated using finite elements (FEM), GMC and RCM are also 
shown in Table [5] 

These results show that the effective two-dimensional bulk modulus is calculated correctly by all the three 
methods. However, the shear moduli calculated for the checkerboard microstructure are quite different from the 
exact result. This exact result also shows that the FEM calculations are the most accurate, followed by GMC and 
then RCM. The values of /Ug ff are also found to most closely approximate the value of G e ff. 

6 Hill's equation 

Hill's equation [7) is an exact relation that is independent of microstructure. This equation is valid for composites 
composed of isotropic components that have the same shear modulus. For a two-dimensional two-component 
composite, this equation can be written as 

' fp H — tF~~p : i (6) 



K eii + G K p + G K b + G 

where / represents a volume fraction, K represents a bulk modulus, and G represents a shear modulus. The sub- 
script 'p' represents a particle property, '6' represents a binder property, and 'eff represents the effective property 
of the composite. 

This relationship is verified using the RVE containing an array of disks occupying 70% of the volume that is 
shown in Figure [5] Table [6] shows the properties of the two components used to compare the predictions of finite 
elements, GMC and RCM with the exact value of bulk modulus predicted by Hill's equation. It should be noted 
that the materials chosen are not quite representative of polymer bonded explosive materials. 

Since the modulus contrast is small, the square array of disks is expected to exhibit nearly isotropic behavior. 
Therefore, the predictions of finite elements, GMC and RCM are expected to be close to the exact values of the 
effective properties of the composite. The numerically calculated values of the effective two-dimensional bulk and 
shear moduli of the composite are shown in Table [7] The percentage difference of the effective bulk modulus from 
the exact value is also shown in the table. 

The effective shear moduli predicted by all the three methods are exact. In case of the effective bulk moduli, 
the RCM predictions are the most accurate followed by GMC and the finite element based calculations. The finite 



Table 4: Comparison of effective moduli for the original and the translated composites. 













Meff 






Meff 






Orig. 


Trans. 


A err (%) 


Orig. 


Trans. 


A OT (%) 


Orig. 


Trans. 


A OT (%) 


FEM 


36.4 


37.8 


-0.8 


10.1 


10 


-3.1 


0.9 


0.9 


22 


RCM 


42.5 


44.5 


6.1 


29.8 


29 


-6.9 


1.3 


1.3 


-292 


GMC 


34.0 


35.1 


-1.3 


3.8 


3.8 


5.3 


0.7 


0.7 


30 
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element based calculations overestimate the effective two-dimensional bulk modulus by around 4.4% while GMC 
underestimates the bulk modulus by around 4.2%. Since the error of estimation of all the three methods is small, 
it is suggested that all three methods are accurate for low contrasts in the shear modulus. However, Hill's equation 
does not appear to be suitable for determining the best numerical method of the three. 

7 Summary and conclusions 

Predictions from the phase interchange identity for the shear modulus are closely approximated by the finite element 
approach (FEM), the recursive cell method (RCM) and the generalized method of cells (GMC) for checkerboard 
composites with low modulus contrast. However, for higher modulus contrasts the FEM approximations of shear 
moduli are the most accurate. The RCM predictions overestimate the shear modulus while GMC underestimates 
the shear modulus. The exact relations for materials that are rigid in shear show that all three numerical techniques 
are inaccurate. The exact relation for this class of materials that is applicable to checkerboard materials shows that 
the FEM calculations are the most accurate while both RCM and GMC perform poorly in comparison. Though the 
predictions of the CLM theorem are quite accurately predicted by all three numerical methods for high modulus 
contrast composites, the FEM results show the least error between the original and the translated effective properties 
while the RCM results show the largest error. The exact relation for isotropic composites with components that 
have the same bulk moduli also shows that the FEM predictions are the most accurate though they are somewhat 
higher than the exact values. However, no such distinction between the three methods can be made using Hill's 
equation. These results reflect previous studies for high modulus contrast, high volume fraction polymer bonded 
explosives using FEM, RCM and RCM and show that exact relations can be used to determine the accuracy of 
numerical methods without performing detailed numerical studies. 
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Appendix 

The components of the two-dimensional stiffness matrix can be computed from two-dimensional plane strain fi- 
nite element analyses. However, the components of the two-dimensional compliance matrix cannot be directly 
determined from two-dimensional plane strain finite element analyses. The reasons for these are discussed in this 
appendix. The approach taken to approximate the two-dimensional compliance matrix is also discussed. 



Table 5: Component properties, exact effective properties and numerically computed effective properties for two- 
component symmetric composite with equal component bulk moduli. 



E v K G 

(xlO 2 ) (xlO 3 ) (xlO 2 ) 

Component 1 25.00 025 2.0 10.0 

Component 2 1.19 0.49 2.0 0.4 

Composite 5.12 0.46 2.0 1.76 





^eff 


Meff 


Diff. 


Meff 


Diff. 




(xlO 2 ) 


(xlO 2 ) 


% 


(xlO 2 ) 


% 


FEM 


20 


1.29 


-26.8 


2.54 


44.4 


GMC 


20 


0.77 


-56.3 


0.77 


-56.3 


RCM 


20 


2.96 


68.0 


4.41 


150.9 
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A.l Two-Dimensional Stiffness Matrix 

The stress-strain relation for an anisotropic linear elastic material is given by 
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For the plane strain assumption, we have, 

£33 = £23 = £31 = 0. (8) 
Therefore, the stress-strain relation can be reduced to 
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The six terms in the apparent two-dimensional stiffness matrix reduce to four is the material is orthotropic, i.e., 



011 




" Cn 


C12 







£11 




C"22 




C12 


C22 







£22 


(10) 


. °12 










Cm 




. e 12 . 





The three constants Cn, C\2 and C22 can be determined by the application of normal displacements in the '1' 
and '2' directions respectively. The constant Cm can be determined using shear displacement boundary conditions 
in a finite element analysis. Hence, it can be seen that the stiffness matrix can be calculated directly from two- 
dimensional plane strain based finite element analyses. This is not true for the compliance matrix. 

A.2 Two-Dimensional Compliance Matrix 

The strain-stress relation for an anisotropic linear elastic material can be written as 
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Table 6: Phase properties used for testing Hill's equation and the exact effective moduli of the composite. 





Vol. 


E 


V 


G 


K 




Frac. 


(xlO 3 ) 




(xlO 3 ) 


(xlO 3 ) 


Disks 


0.7 


3.00 


0.25 


1.20 


2.40 


Binder 


0.3 


3.58 


0.49 


1.20 


60.00 


Composite 


1.0 


3.22 


0.34 


1.20 


3.82 
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Table 7: Numerically computed effective 



properties for a square array of disks with equal component shear moduli. 





Kcs 


% Diff. 




Meff 




(xlO 3 ) 




(xlO 3 ) 


(xlO 3 ) 


FEM 


3.98 


4.4 


1.20 


1.20 


RCM 


3.92 


2.7 


1.20 


1.20 


GMC 


3.66 


-4.2 


1.20 


1.20 



The relationship between the stiffness matrix and the compliance matrix is 
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C66 



(12) 



or, 



S = C 



(13) 



It is obvious from the above equation that the apparent two-dimensional compliance matrix is not equal to the 
inverse of the apparent two-dimensional stiffness matrix, i.e., 



511 S12 S\q 

512 S22 S26 

5*16 S26 Sqq 



Cn C\2 Ci6 
C\2 C22 C26 
Cl6 C16 C66 



(14) 



Hence, we cannot determine the two-dimensional compliance matrix if we only know the two-dimensional stiffness 
matrix. 

Let us again examine the effect of the plane-strain assumption on the stress-strain relation. We then have 



(15) 





en 






S12 


5*13 


Sw 




Oil 




£22 




•512 


S22 


^23 


S26 




C22 









Sl3 


S23 


5*33 


S36 




C"33 




. £ 12 . 




Sie 


5*26 


^36 


See 




_ (T12 


For orthotropic materials, this relation simplifies to 












en 




' S u 


Sl2 


5*13 







o-n 




£22 




S12 


S22 


£23 







0"22 









S13 


S23 


S33 







C"33 




. £ 12 . 













See 




_ (T12 



(16) 



This equation shows that we need to know the stress 033 to determine the terms of the compliance matrix and hence 
three-dimensional analyses are necessary. If we assume plane stress, we can determine the terms of the matrix S 
directly. However, the apparent two-dimensional compliance matrix for plane stress is not equal to that for plane 
strain and hence we cannot apply this method for our purposes. This is why the plane strain compliance matrix 
cannot be determined using two-dimensional finite element analyses only. 
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A.3 Approximation of Compliance Matrix 



The two-dimensional compliance matrix can be determined approximately for materials with square symmetry by 
assuming that 513, S23 and S33 are known. Let, 

v z 1 

Sl3 = 5*23 = — — S33 = — (17) 
^3 ^3 

where, v 3 is the Poisson's ratio in the out-of-plane direction and E% is the Young's ratio in that direction. Then, for 
a material with square symmetry, 

en 

£22 



Inverting the relation, we have, 

CTll 

C22 
033 

0"12 



S\2 

"3 

E 3 





612 
Sii 

"3 

E 3 





^3 

Ez 

I 3 









Sm 









022 




0"33 




0"12 



(18) 



Cn 


C12 


C13 







en 


C12 


Cn 


C23 







£22 


C13 


C23 


C33 



















Cm 




£12 



(19) 



whrrr C\, — EzSxx-vl 

wnere, Uu- EsSf^qs^-EsS^^s^ , 
C12 = 



B35f 1 -2^5ii- J B3S'f 2 +2y 2 Si2 ' 



Note that it is not necessary to know C13, C23 and C33 to determine Sn and Si 2 . 
We can write the above relations between Cn, C12 and Sn, Si 2 in the form 



E3SI " + 2^) Sn - (e 3 S 2 2 - MS12 ~^-) =0, (20) 

E 3 S 2 12 ~(^- + 2i%\ S12 - (e 3 S 2 u - 2u 2 3 S u + = 0. (21) 



In simplified form, 



12 / V W2 



AiS u + Si5n + Ci = 0, (22) 
A 2 Si 2 + B 2 S 12 + C 2 = 0. (23) 



We can solve these quadratic equations to get expressions for Sn and Si 2 as 

511 =^±f^, (24) 

512 = - B -^ MU . (25) 

Knowing Cn, C12, £3 and v 3 these two equations can be solved iteratively to determine Sn and £12. The values 
of Cn and C12 can be determined using the procedure outlined at the beginning of this section. It remains to be 
discussed how £3 and v 3 are to be determined. 
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A.4 Determination of E 3 and i/ 3 



Two methods can be used to determine the values of E3 and 1/3 for our calculations. The first method is to assume 
that the rule of mixtures is accurate enough to determine the effective properties in the '3' direction. Thus, if the 
volume fraction of the first component is f\ and that of the second component is f 2 , we have, 

E 3 =fiE 1 + f 2 E 2 , (26) 

^3 =/l^l+/2^2- (27) 

where Ei and Vi are the Young's modulus and the Poisson's ratio of the ith component. 

The other option is to use the values of S13, S23 and S33 obtained from GMC since these are also quite accurate 
for the out of plane direction. Thus, we have, 

E3 = 5<Jmc> (28) 

^3 =-^3 MC ^3- (29) 

This is the procedure we have use to determine the effective compliance matrices discussed in this paper. 
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Figure 5: RVE for a square array of disks. 
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